function [x,y,z,fitresults,invarx,invary,invarz,graymask,whitemask]=marmplots(xin,yin,zin,maskselection,figurenum,subplotpos)
% function [x,y,z,fitobject,invarx,invary,invarz,graymask,whitemask]=marmplots(xin,yin,zin,maskselection,figurenum,subplotpos)

%% setup outputs
x = [];
y = [];
z = [];
fitresults = [];
invarx = [];
invary = [];
invarz = [];
graymask = [];
whitemask = [];

%% set RECOMMENDED input system parameters
if nargin<6
    subplotpos = 1:100;
end
if nargin<5
    figurenum = 1;
end

%% load masks
marmmasksetter;

%% gather x data
marmxdatasetter;

%% gather y data
marmydatasetter;

%% gather z data
marmzdatasetter;

%% set x axis label
marmxaxislabelsetter;

%% set y axis label
marmyaxislabelsetter;

%% set z axis label
marmzaxislabelsetter;

%% process data and plot results
marmdataprocessorandplotter;


